Vibration-induced correction to the current through a single molecule 
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We provide analytical results for the perturbative correction to the current-voltage relation 
through a vibrating molecule for weak electron-phonon coupling. The nonlinear conductance ex- 
hibits a steplike feature at eV = fiwo, where Wq is the vibration frequency. We establish criteria 
for the sign change of the step in the conductance (up or down). This transition turns out to be 
nonuniversal and is governed by essentially all system parameters. 
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Electronic transport through single molecules has at- 
tracted much attention lately; for recent reviews, see 
Rcfs. UIllSHS- Besides the technological promises 
raised by molecular electronics, this field also poses in- 
teresting questions to theory. In this Brief Report, 
we will revisit the problem of how the coupling to 
a vibrational mode ("phonon") of frequency loq and 
electron-phonon coupling strength g affects the cur- 
rent through a single molecule. The resulting fea- 
tures in the I-V characteristics are often referred to 
as inelastic electron tunneling spectroscopy^ Theories 
for various aspects of this rich and diverse problem 
have been proposed over the past few years by a 
large number of authors ^ 7 -^ ' 10 i 11 ' 12 i 13 primarily mo- 
tivated by groundbreaking experiments demonstrating 
the influence of vibrational degrees of freedom in single- 
molecule transport. To mention just a few key exper- 
iments, single-molecule transport has been studied us- 
ing various organic molecules^ fullerenes , 1 5 1 16 1 1 7 carbon 
nanotubeS ) 18 i 19 or single hydrogen molecules between Pt 
leads' Phonon-assisted processes were shown to imply a 
step in the I-V characteristics once the dc bias V reaches 
the threshold value huo/e for excitation of a phonon 
mode. Related features can sometimes be seen at inte- 
ger multiples of this value. Vibrational effects on single- 
molecule transport have recently been reviewed^ includ- 
ing a discussion of the validity regime for our Hamiltonian 
below. 

It is remarkable that the experimentally observed 
vibration-induced step features in the differential con- 
ductance can either decrease^- or increase the conduc- 
tance through the molecule near eV = fc^o' 14115 ' 171 " 18 
This corresponds to dips versus peaks in the second 
derivative, d 2 I/dV 2 . Such features have also been re- 
ported theoretically 5 ' 21 ^ 2 Recent theoretical work on 
this questional!! argues that the transparency T of the 
single-molecule junction basically determines the step di- 
rection. The critical value was reported to be T = 1/2, 
with a step down (up) in the differential conductance 
at eV = hu for T > 1/2 (T < 1/2). This conclu- 
sion seems roughly consistent with existing experimental 
data: For the measurements^ T close to unity was 
reported, while typically T < 0.1 in the other experi- 



ments. However, given the many parameters present in 
even the simplest Hamiltonian, one may question why 
the crossover should be universal in the sense that it is 
only determined by the transparency T = 1/2. Here, 
we reexamine the question of current increase or de- 
crease at the phonon excitation threshold eV = Hluq. 
We derive and discuss analytical results for the current 
correction SI perturbative in the electron-phonon cou- 
pling strength g. Experimental values for g are often 
very small, justifying a truncation of the perturbation 
series already at the lowest nontrivial order. For that 
reason, our expressions below are expected to provide 
useful estimates for many experiments. However, we 
will not attempt a detailed description of specific exper- 
iments, but instead aim at an analytical understanding 
of vibrational features in the I-V characteristics under 
a simple yet realistic model. In our opinion, a thor- 
ough understanding of the lowest-order feature is worth- 
while, given the complexity of the physical processes in- 
volved. Most available results are obtained from lengthy 
numerical calculations and do not easily yield general 
insights. In addition, many calculations were based on 
essentially uncontrolled approximations, rendering their 
predictive power questionable. Moreover, some published 
theoretical work on this subject have used approximate 
schemes that are in conflict with the basic requirement 
of current conservation. Current conservation is auto- 
matically fulfilled under self-consistent approximations, 
but is generally violated otherwise (unless particle-hole 
symmetry is present), limiting the practical usefulness of 
many approximations to special parameter sets^ The 
lowest-order correction, however, can be evaluated ex- 
actly and, therefore, does not suffer from any such limi- 
tations. Although aspects of the lowest-order correction 
51 have been studied before ^ 21 ' 22 a complete and analyt- 
ical discussion was not given so far. We show that besides 
the step feature caused by inelastic processes, quasielastic 
electron-phonon interactions are responsible for another 
singular term that logarithmically diverges at eV = HajQ. 



We study the model of just one relevant molecular level 
( "dot" ) , the so-called local Holstcin model, also employed 
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by previous work^ 

H = (e + gQ)d ] 'd + tuu b' i b + 
+ (^c ka + H.c.) 
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which neglects the Coulomb interaction U and is formu- 
lated for spinless dot fermion operator d (we set h = 1). 
The effect of the lowest-order correction in U, consistent 
with our g 2 calculation below, is anyway trivial and can 
be absorbed by a renormalization of the chemical poten- 
tial. We take the standard wide-band limit for the leads, 
which is justified if the lead density of states v$ does not 
vary significantly in energy on the relevant scales. The 
lead modes are occupied according to Fermi functions 
/L/fl(e) = /(e - Vl/r), where \i L - fi R = eV defines 
the bias voltage V. We introduce the hybridizations as 
Flm. = nv \j L / R \ 2 , and define r = T L + T R . Finally, 
the boson operator b describes an Einstein phonon mode 
(vibration mode) of frequency lu , with linear coupling of 
strength g between the displacement operator Q = b + b^ 
and the dot occupation operator d'd. The electrical cur- 
rent through the dot can be computed from the retarded 
dot Green's function (GF) evaluated in the presence of 
the leads and the phonon, G r (u>), according to the well- 
known expression^! 
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The non-interacting (g = 0) Keldysh GF describing the 
out-of-equilibrium dot coupled to the leads is 
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Note that we use the unrotated Keldysh notation, where 
the connection to retarded or advanced GFs G r l a and 
the lesser GF G < is established by 



G- 



G-+ 
G++ 



G r + G< G< 
G r -G a + G< -G a + G< 



such that Gq(u>) — (ui— eo+iT)^ 1 . Since G T obeys its own 
Dyson equation, we only need to compute the retarded 
self-energy Y7(lS) to order g 2 , resulting in 



G r {u) = Gl^) + Gl{u)Tr{u)Gl^), 



(4) 



where the second term defines the correction 51 in 
Eq. @ • We will focus on the most interesting T = limit 
from now on, where the Fermi function is f(u>) = 0(— u>). 
Defining 
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eo, V a =L/R=± = M ± eV/2, (5) 



the first (g = 0) term yields 
e 4T L T R 



Io(V) = -- 
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[tan" 1 (ml /r) - tan- 1 (jl R /T)] . (6) 



The V — > transparency of the junction, T 
(h/e 2 )dI/dV, follows as 



4r L r, 
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Note that eo an d the mean chemical potential always 
appear in the combined scale /2. 

Let us now analyze the lowest-order correction 51 to 
the current. It arises from the retarded self-energy S r (w) 
evaluated to order g 2 due to phonon processes, entering 
Eqs. (TJ| and @. There arc two contributions coming 
from the real (imaginary) parts Y7 R (EJ), corresponding 
to quasielastic (inelastic) processes, 
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The self-energy is readily computed on the perturbative 
level, where two diagrams are present in order g 2 . The 
"tadpole" diagram docs not carry frequency dependence 
and can be absorbed by a renormalization of p,. We there- 
fore keep only the standard "rainbow" diagram, which 
gives the retarded self-energy E r (tj) = E (w) + £ h (cj) 
from 



: (cj) = Tig' 
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D - ± (fi)G ± (w-fi). (10) 



Here, Dq(u>) is the 5 = GF of the displacement operator 
Q, which for T = is given by 
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Using a Wiener-Hopf decomposition of Gq~ (lu) into the 
parts analytic in the upper or lower complex frequency 
plane, one arrives at the result (cf. also Ref. 0), 
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The computation of <57 is then reduced to a single fre- 
quency integration. We see that due to the phonon mode, 
the retarded electron self-energy contains directly the 
Fermi functions shifted by ±luq. Hence, a singular (step) 
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dependence of its imaginary part on the energy results, 
which must be accompanied by a logarithmic singularity 
in the real part due to analytic properties. These sin- 
gularities have been discovered first by Engelsberg and 
Schrieffer in their study of bulk Einstein phonons^ With 
the above self-energies, it is easy to check that the re- 
quirement for current conservation^ 



duj [G < (cj)E > (w) - G > (w)S < (cj)] = 0, 



is fulfilled (also at finite T) to the required g 2 order. 
Let us first discuss the inelastic part, 5I mc \. Using the 
auxiliary relation for an arbitrary function F{u>, V), 



d f^ L e f fJ ' L 

— / cLjF(uj, V) = -J2 F(fi ai V)+ / du 



dF(u, V) 
dV ' 
(14) 

some algebra gives the g 2 inelastic correction to the T = 
nonlinear conductance for arbitrary parameters in closed 
form, 



dSI h 
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We focus now on the singular contribution at eV = Huo, 
which is best illustrated by analyzing d 2 5I/dV 2 . Singular 
terms come from the derivative acting on the Heaviside 
function in Eq. (fT5|). and produce a delta peak, 
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with the dimcnsionlcss coefficient 

1 - [(/2/r) 2 - K/2r) 2 ] 2 + 2lu (T l - T R )fl/T 3 
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(17) 

which is the main result of this Brief Report. We will 
discuss this result below, but first turn to the quasielas- 
tic contribution 5I qc i due to the real part (|12p of the 
self-energy. Using Eq. (|14j) , we find again a singular con- 
tribution in the differential conductance at V = hwo/e. 
We obtain the analytical result, valid for V ~ hu>o/e, 
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with the dimensionlcss coefficient 
c q ci = 
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[1 + (p/r + ^ /2r) 2 ] 2 [i + (fi/r - Wo /2r) 2 ] 2 ' 

(19) 



It is worth mentioning that c qc i = for large r > wo 
at the symmetric point Tl = I\r. Away from this 
limit, however, the logarithmic term in Eq. (|18| will be 
present. All other contributions to <5/ qe i beyond Eq. (TT8"]) 
are smooth and featureless at eV « hujQ, and do not af- 
fect the characteristic feature in d 2 I/dV 2 , whereas the 
singular contribution (|18[) is logarithmically divergent at 
the threshold. Note that this logarithmic divergence due 
to quasi-elastic processes creates a symmetric dip or peak 
(depending on the sign of c qe i) in the differential conduc- 
tance at eV = hujQ, while the inelastic contributions are 
responsible for a step feature. In the full d 2 I/dV 2 curve, 
this translates to asymmetric dips or peaks. The relative 
importance of inelastic versus quasi-elastic contributions 
can be judged from the ratio Ci no i/c qo i. For the symmet- 
ric case, = Tr = r/2, a simple result follows from 
Eqs. (HZJ and (HI), 
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For small loq/T or large /2, we have r> 1 and the inelastic 
channel always dominates, while for large luq, quasielas- 
tic processes can be more important. The perturbative 
results (fT8| and (fT6|) obviously break down close to the 
threshold voltage. At T — and in the absence of an 
external bath, the only way to account for the finite life- 
time of the phonon, and hence the smearing of the step 
and/or peak features, is to take into account the elec- 
tronic polarization in the phonon GF. The retarded po- 
larization function x r ( w ) W1 H then result in a damping 
7 ~ — g 2 Xi( UJ o) of the phonon mode. We obtain after 
some algebra the nonequilibrium electronic polarization 
function in analytical form. In the particle-hole symmet- 
ric case, this result simplifies to 



x r H = - 



i 



7r lu(lo + 2iT) 



In 



V r 2 + u 2 /4y v ' 



This implies the estimate 7 ~ g 2 ujo/irT 2 in the limit 
of a soft phonon loq T, and 7 ~ g 2 T /us 2 , for a hard 
phonon luq 3> T. However, phonon damping is, in fact, a 
higher-order effect in the electron-phonon coupling, and 
to consistently account for the finite damping 7 while 
keeping current conservation intact remains a theoretical 
challenge. Other effects of higher-order diagrams include 
the proliferation of steps and/or peaks at multiples of 
the phonon frequency ujq. Indeed, it is easy to see that 
the g 2 "th-order rainbow diagram in the electronic self- 
energy produces a step feature in the differential conduc- 
tance at the voltage V — ntuvo/e. The appearance of 
such step features at multiples of ujq is closely related to 
the strong-coupling picture obtained through a polaron 
transformation^ However, when going beyond the low- 
est order in g, vertex corrections are also expected to 
be important. Unfortunately, the proper treatment of 
such nonequilibrium many-body effects (respecting the 
requirements posed by current conservation) remains a 
challenging task and is beyond the scope of this work. For 
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weak electron-phonon coupling, which appears to be ap- 
propriate for many experiments, none of these fine details 
matter in any case, and the g 2 calculation is sufficient. 
The d 2 I/dV 2 feature is then determined by Eqs. (Til)]) 
and (|18p , where the damping 7 in the phonon mode acts 
to broaden the delta function in Eq. (|16p within a phe- 
nomenological description. 

The above results allow us to clarify the question of 
peak vs dip in the second derivative, d 2 I/dV 2 , which 
arises due to the singular inelastic correction (|16p . For 
= Tr and p, = 0, where the transparency (|7|) is ideal, 
T = 1, we observe from Eq. (fT?)) that for u > 2r, in- 
stead of the expected dip (ci no i > 0), one actually ob- 
serves a peak. For /i ^ 0, once \p,\ > \JY 2 + Wg/4, one 
finds a peak. Note that for p, = ±r, the transparency 
([7]) is precisely 1/2, thereby allowing us to rationalize 
why previous numerical studies for related models^^ re- 
ported a T = 1/2 criterion for the transition from peak 
to dip. This value correctly describes the transition in 
the limit of a soft phonon, luq <C T, and assuming sym- 
metric contacts, Tl = Tr. The value T = 1/2 was, in 
fact, established precisely in this parameter regime^^ 
Our analytical result (jTTJ) , shows, however, that, in gen- 



eral, the transition is nonunivcrsal and determined by 
all parameters. For example, it can be achieved cither 
by tuning T — where the precise transition value de- 
pends also on luq and the asymmetry Tl — Tr, and is 
only approximately given by T = 1/2 — or by changing 
other parameters, such as Tr - or the ratio between 
phonon frequency and hybridization, ujq/T. The nonuni- 
versality of the step is also implicit in the early work 
on phonon-assisted tunneling through a resonant level 
by Glazman and Shekhter^l On top of this peak or dip 
structure due to inelastic processes, the quasielastic con- 
tribution causes a singular response near the threshold 
value eV = Hluq. This logarithmic correction to the dif- 
ferential conductance implies an asymmetric line shape 
in d 2 I/dV 2 as discussed above. Such asymmetries have 
frequently been reported experimentally and in numer- 
ical calculations,— and they are a direct consequence of 
the Engclsbcrg-Schricffer singularity. 
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